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Abstract 

The objective of the present study was to estimate the level of inbreeding coefficient during 
inbreeding of the pedigree of random mating and selected populations of two distinct races of 
mulberry silkworm, Bombyx mori (Lepidoptera: Bombycidae), in the silkworm germplasm. The 
six generation data of the two races, namely multivoltine Pure Mysore and bivoltine NB4D2, were 
studied for inbreeding depression coefficient using the residual maximum likelihood method, 
utilizing two statistical models by analyzing six quantitative traits, namely, larval weight, cocoon 
weight, shell weight, shell ratio, pupation rate, and filament length. The results of the present 
experiment demonstrated that the inbreeding coefficient was significant in Model 1 for most of 
the economic traits in the random mating populations of both the races compared to those of 
selected populations. These results suggest that during stock maintenance, application of rigid 
selection for increased numbers of generations helps to retain original characteristics of the pure 
races while reducing the deleterious effects of inbreeding. The significance of inbreeding 
coefficient is discussed with reference to the inbreeding of silk moths in the silkworm 
germplasm. 
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Introduction 



The mulberry silkworm, Bombyx mori L 
(Lepidoptera: Bombycidae), is a lepidopteron 
insect known for the production of silk 
("queen of natural fibers"), and extensive 
genetic studies have been carried out utilizing 
this insect (Tazima 1978). Sericulture, an 
important source of revenue generation in 
many countries, is based on the maintenance 
of pure races belonging to different voltinistic 
groups in germplasm centers. Sib-mating of 
the progenies derived from a single brood is 
preferred, and original traits of the silkworm 
races are protected to a great extent. The 
homozygous races maintained in the 
germplasm centers are subsequently crossed 
to produce hybrids for commercial 
exploitation. As a result, silkworm breeders 
retain the original characters of the races 
through selection pressure to avoid inbreeding 
depression. 

The inbreeding depression for a trait is a linear 
function of the inbreeding coefficient (Lynch 
and Walsh 1998). One of the measures of 
inbreeding is the inbreeding depression 
coefficient, which corresponds to half of the 
relation between parents. Utilization of the 
linear-regression method to estimate 
inbreeding depression for production traits in 
livestock breeding is well documented 
(Allaire and Henderson 1965; Henderson 
1984; Hudson and Van Vleck 1984; Wall et 
al. 2005). In sand cricket, Gryllus firmus, it is 
reported that inbreeding depression can cause 
a substantial decrease in trait values, and life 
history traits show significantly higher 
inbreeding depression than morphological 
traits (Roff 1998). In mulberry silkworm, 
there are a few reports on the inbreeding 
coefficient for economic traits of cocoons 
(Nematollahian et al. 2010; Azizian et al. 
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2011). The problem of applying inbreeding to 
silkworm population has been discussed by 
Petkov et al. (1999). However, there is no 
information available in regard to the effects 
of selection on inbreeding depression 
coefficient for various traits in germplasm 
centers. Thus, the purpose of this research was 
to study (1) the level of inbreeding depression 
coefficient for six quantitative traits in two 
races inbred for six generations and (2) the 
differential expression of inbreeding 
coefficient for two different models. 

Materials and Methods 

In order to estimate the level of inbreeding 
and investigate the occurrence of inbreeding 
depression through inbreeding coefficient in 
the random mating (non-selected) and selected 
populations, a pilot study was undertaken by 
drawing two distinct races of the mulberry 
silkworm, namely multivoltine Pure Mysore 
and bivoltine NB4D2, from the germplasm 
bank of the Department of Studies in 
Sericulture Science, Mysore, India. The 
parental population of the two races were 
grouped into random mating and selected 
populations. In the random mating 
populations, selection was not applied for the 
six traits under study, whereas for the selected 
populations, selection pressure was applied at 
every stage of the life cycle of silk moth (egg, 
larva, pupa, and adult) comprising six 
quantitative traits (larval weight, cocoon 
weight, shell weight, shell percentage, 
pupation rate, and filament length), which 
were selected during rearing of silkworm for 
six consecutive generations. The two 
populations were maintained in replicates of 
three for statistical analysis. The results were 
computed according to the equation described 
by Falconer (1989), and the single trait linear 
equation was applied as detailed below. 
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Rate of Inbreeding per generation: F t = AF + 
(l-AF)Ft-i 

Where: 

F t = inbreeding at generation t, F t _i = 
inbreeding at generation t-1, AF = inbreeding 
coefficient. 

The next step is the estimation of the AF from 
the above formula. The rate of inbreeding 
(AF) was estimated as the difference between 
the inbreeding of the individual (F t ) and the 
average inbreeding of the parents (F t _i) 
divided by (1- F t -i). (Falconer and Mackay 
1996). 

AF = F t - F t -i / 1 - F t -i 

When the inbreeding coefficient is expressed 
in terms of AF, the above formula is valid for 
any breeding system and is not restricted to 
the idealized population. 
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quantitative trait was included as one response 
variable and fixed effect, which is a covariate 
and includes interaction between races, 
replicates generation, and is a constant value. 
In Model 2, in addition to the above, 
generation squared (G 2 ) was added to the 
fixed effects. 

Model 1 = constant + races + races* replicates 
+ AF + generations 

Model 2 = constant + races + races* replicates 
+ AF + generations + generations 

The formula is as follows: 



Vr 
U 




(1 + F A ) 







Where F x is the coefficient of inbreeding of 
individual x, n is the number of connecting 
links between the two parents of x through 
common ancestors, and Fa is the coefficient of 
inbreeding of the common ancestor. 



For repeated full sib-mating, the inbreeding 
coefficient settles down to a constant value of 
0.191/generation. Hence, AF = 
0.191/generation (Falconer 1989). Based on 
the above equation, the data is subjected to 
two linear mixed-models, using the residual 
maximum likelihood method directive in 
GENSTAT 9 (2006). In Model 1, each 



Results and Discussion 

The mean values and inbreeding coefficient 
(AF values) for the two races, both in the 
random mating and selected populations of 
silkworm, are listed in Table 1 , and the results 
are herein discussed. 



Table I . Inbreeding coefficient for six quantitative traits in the random mating and selected populations of the two races of the 
mulberry silkworm, Bombyx mori. 



Race 


Trait 


Random populations 


Selected populations 


Mean 


Model 1 (AF) 


Model 2 (AF) 


Mean 


Model 1 (AF) 


Model 2 (AF) 




Larval weight (g) 


42.870 ±0.5 11 


-41.500 ± 18.525 


-152.800 ± 102.410 


44.960 ± 0.345 


-36.090 ± 12.488 


-54.950 ± 138.066 


Bivoltine NB4D2 


Cocoon weight (g) 


1.839 ±0.012 


-0.796 ± 0.423 


4.965 = 4.365 


1.844 ±0.011 


-1.070 ±0.384 


3.016 = 4.079 


Shell weight (g) 


0.324 ± 0.004 


-0.354 ±0.148 


0.654 ±0.616 


0.345 ± 0.004 


-0.433 ±0.150 


2.116 ± 2.481 


Shell ratio (%) 


17.970 ±0.1 44 


-12.260 ±5.224 


-64.140 ±69.973 


18.780 ±0.1 13 


-11. 900 ±4.099 


57.640 ± 60.625 


Filament length (m) 


837.800 ± 11.640 


-767.900 ±421.550 


-5931 ±4417.300 


872.800 ± 11.310 


-1136*409.700 


3953 ± 4285.9 


Pupation rate (%) 


91.720 ±0.180 


-7.639 ±6.533 


64.000 ±53.401 


93.390 ±0.175 


-29.310 ±6.328 


-105.000 ± 166.480 


Multivoltine PM 


Larval weight (g) 


21.190 ±0.234 


-23.270 x 8.475 


-9.444 ±93.679 


23.220 x 0.223 


-20.690 X 8.078 


-50.310X83.402 


Cocoon weight (g) 
Shell weight(g) 


1.008 -T- 0.005 
0.120 ^ 0.001 


-0.258 -r 0.187 
-0.028 ±0.019 


-5.207 -T- 1.498 
-0.891 ± 0.213 


1.094 -T- 0.003 
0.141 ±0.002 


-0.323 -0.127 
-0.102 = 0.043 


-0.366-0.423 
-0.306 ±0.475 


Shell ratio (%) 


11. 930 ±0.031 


1.666 ±1.137 


-26.590 ± 9.549 


12.860 ±0.081 


-6.631 ±2.919 


-42.580 ± 50.572 


Filament length (m) 


324.300 ±1.920 


-73.277 ± 69.595 


-1036 ± 709.400 


341.800 ± 1.940 


-141.1 ±70.350 


-336.000 ± 776.27 


Pupation rate (%) 


93.000 + 0.128 


-13.010 ±4.641 


-77.650 ±5 1.332 


95.50 ±0.156 


-12.660 ±5.666 


-21.180 ±61.911 



Model I = Constant + Races + Races * Replicates + AF + Generations 

Model 2 = Constant + Races + Races * Replicates + AF + Generations + Generations squared 
AF = Inbreeding coefficient 
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Random mating populations 

In Model 1, all six analyzed traits in the 
random populations exhibited negative 
significant AF values where generation (G) 
was fitted as one of the covariates. In the 
bivoltine NB4D2 race, the least AF value being 
observed for the trait shell weight (- 0.354) 
and the highest for filament length (- 767.9) 
where as in multivoltine Pure Mysore race the 
AF value ranges from a lowest of -0.028 for 
shell weight to a highest of -73.277 for the 
trait filament length. 

In Model 2, where G 2 was fitted as one of the 
covariates, two of the six traits, namely larval 
weight and filament length, analayzed in the 
NB4D2 race showed significant negative AF 
values (-152.800 ± 102.410 for larval weight 
and 

-5931 ± 4417.3 for filament length), and a 
non-significant negative AF value was 
observed only for one trait, namely shell ratio 
(-64.140 ± 69.973). Contrary to the above, the 
other three traits, namely cocoon weight 
(4.965 ± 4.365), shell weight (0.654 ± 0.616), 
and pupation rate (64.000 ± 53.401), exhibited 
respective positive significant AF values. In 
the case of the multivoltine Pure Mysore race, 
five of the six traits fell into the negative 
significant group, while larval weight fell ino 
the negative non-significant group (-9.444 ± 
93.679). The lowest AF values appeared for 
cocoon weight and shell weight in both the 
races. In the progenies of the random mating 
populations, the low inbreeding depression 
coefficient values (e.g., cocoon weight and 
shell weight) in both the models clearly 
indicate low dominance variance for these 
traits, and vice versa is true for progenies of 
selected populations (e.g., filament length). 
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Selected populations 

In the selected populations of the two races 
studied in Model 1, all six traits showed 
negative significant AF values. In the 
bivoltine NB 4 D 2 race, the lowest and highest 
AF values were -0.433 for shell weight and - 
1136 for filament length, respectively. 
Similarly in the Pure Mysore race, the lowest 
AF value (-0.102) was recorded for shell 
weight, and the highest AF value (-141.1) was 
recorded for filament length. The traits cocoon 
weight and shell weight exhibited low AF 
values in both the races under study. 

In Model 2, all six traits selected for analysis 
recorded non-significant AF values for both 
the races. In the NB 4 D 2 race, larval weight (- 
54.950 ± 138.066) and pupation rate (- 
105.000 ± 166.480) showed negative non 
significant AF values, whereas the remaining 
four traits (cocoon weight, shell weight, shell 
ratio, and filament length) exhibited 
respective positive significant AF values of 
3.016 ± 4.079, 2.116 ± 2.481, 57.640 ± 
60.625, and 3953 ± 4285.9. The lowest value 
noticed was for shell weight (2.116), and the 
highest value noticed was for filament length 
(3953). In the Pure Mysore race, all the traits 
showed negative non-significant AF values. 
The lowest value was recorded for shell 
weight (-0.306) and followed by cocoon 
weight (-0.366). The highest value was for 
filament length (-336.0). Cocoon weight and 
shell weight had the lowest AF values in both 
the races. It is quite interesting to note that 
both the lowest (shell weight) and highest 
(filament length) inbreeding depression 
coefficient values were noticed for the same 
traits in random mating and selected 
populations of the two races. From the results 
it is apparent that when Model 1 was applied, 
irrespective of the type of the population and 
the involvement of selection pressure, all the 
traits showed significant inbreeding 
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depression. However, when Model 2 was 
applied, the results clearly demonstrated that 
in random mating populations most of the 
traits showed significant inbreeding 
depression values, except shell ratio in 
bivoltine NB4D2 race and larval weight in 
Pure Mysore race. In selected populations for 
the same model, all the traits being studied 
revealed non-significant inbreeding 
depression values. The absence of inbreeding 
depression in selected populations for all the 
traits studied is an indicative of the dominant 
gene effect. 

The analysis of the results by applying the two 
models clearly indicates that during stock 
maintenance in germplasm centers, 
application of rigid selection for more 
numbers of generations helps in the retention 
of original characteristics of the pure races 
and also reduces the deleterious effects of 
inbreeding. Thus, the net effort of inbreeding 
in the selection program will depend on the 
magnitude of the selection response relative to 
the possible depression and rate of 
accumulation of inbreeding, as was reported 
by Swanepoel et al. (2007) in sheep. A 
comparative account of the values of AF 
between the two populations being studied 
clearly revealed low AF values for cocoon 
weight and shell weight, two traits that are 
important to silk production. Therefore, if 
variability is not introduced through selection, 
there are possibilities of increases in 
inbreeding coefficient values, and the 
deterioration of these traits are expected. 

In conclusion, it is suggested that the breeders 
of germplasm take notice of AF values. 
Prolonged maintenance of stocks (G 2 ) with 
rigid selection, rather than short duration (only 
G) of selection, will be useful in race 
maintenance programs to minimize the 
inbreeding depression. 
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